Radiative coupling of quantum dots in photonic crystal structures 
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We derive a general formalism to model the polariton states resulting from the radiation-matter 
interaction between an arbitrary number of excitonic transitions in semiconductor quantum dots 
and photon modes in a photonic crystal structure in which the quantum dots are embedded. The 
Maxwell equations, including the linear nonlocal susceptibility of the exciton transitions in the 
quantum dots, are cast into an eigenvalue problem, where the matrix elements are expressed as 
functions of the microscopic parameters of the semiconductor system. We compute realistic photon 
modes using Bloch-mode expansion, while the simplest assumption of a delta-like excitonic wave- 
function for the quantum dots is made. As examples, we apply the formalism to the case of one 
or two quantum dots embedded in a photonic crystal with a line defect - either an Ln cavity or 
a Wl waveguide. For the single dot case, we reproduce known analytical results, while for the 
two dot case, we study the radiative excitation transfer mechanism and characterize its strength, 
the dependence on the detuning between quantum dot and photon modes, and the dependence on 
inter-dot distance. We find in particular that, for typical InGaAs quantum dots in a GaAs photonic 
crystal structure, the inter-dot radiative coupling strength can reach 100 /xeV in a short cavity, and 
its decay with distance in longer cavities and waveguides is determined by the group velocity of the 
exchanged photons and their radiative lifetime. We finally show that, for an Ln cavity of increasing 
length, the radiative excitation transfer mechanism is subject to a crossover from a regime where 
a single photon mode dominates the transfer process, to a multi-mode regime. For the specific 
system under study, this crossover occurs around n = 150. The numerical procedure developed in 
the present work can in principle be applied to any structure whose photon modes can be computed 
with reliable accuracy, and in addition naturally allows for disorder effects in the photonic crystal 
structure to be taken into account. 

PACS numbers: 78.67.Hc, 71.36.-l-c, 42.50.Ct, 42.70.Qs, 03.67.-a 
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I. INTRODUCTION 



In the last two decades, the design and implementation 
of devices for quantum information processing has been 
a major goal of condensed matter physics. An essen- 
tial requirement of the quantum information paradigm 
is the possibility for two qubits to interact coherently in 
a controlled fashion, in order to achieve controlled gate 
operations. This must in principle be possible for each 
arbitrarily chosen pair of qubits in the system. Most of 
the technologies, however, employ qubits which are at all 
times spatially separated and do not interact directlji"— . 
The interaction can then be achieved by means of a quan- 
tum bus, namely a spatially extended degree of freedom 
interacting with all localized qubits. In a more gen- 
eral picture, these spatially extended degrees of freedom 
might even form a quantum network connecting distant 
quantum information systems A quantum bus can be 
of several kinds - two common examples being phonons 
in chains of trapped ions^ and microwave photons in su- 
perconducting circuitsjS'^. Photons are the most natural 
choice in a solid-state system, given their low decoher- 
ence rate, high velocity, and the recent advances in the 
on-chip photonic technology, especially in the photonic 
crystal (PHC) domain^ii^. There, extremely high-Q cavi- 
ties with modal volumes of the order of (A/n)^ have been 
fabricated both in silico n^°d^ and in GaAsi^, as well 
as waveguides allowing for low-loss, long-range photon 
transfer with a controllable group velocity^iS The advance 



in PHC technology opened the way to several exper- 
imental breakthroughs, including low power all-optical 
switchingfi^^ and the dynamic control of the strong 
coupling between two distant cavities^ - highlighting 
the extreme level of control of light which is currently 
achievable. 

Semiconductor quantum dots (QDs) have long been 
considered as viable qubit candidates,— as they natu- 
rally fulfill the criteria of scalability and integrability re- 
quired in a quantum information technology. Facing the 
remarkable advance made in the system of spin qubits 
in lateral QDs^. - where electron spins are controlled 
electronically with ohmic contacts - optical excitations 
in self-organized QDs have only recently caught up in 
the race towards controlled quantum operations. On 
one hand, in fact, full singie-qubit optical control has 
been successfully demonstrated J^"— On the other, in- 
tegrating QDs in photonic structures has made signifi- 
cant progress, and both single-dot Purcell enhancement 
in cavities^^"'^^ and waveguides^"—, and strong coupling 
to a cavity mode^^"~ have been demonstrated. Single- 
dot coupling to light modes is in itself important for 
practical applications, as suggested by the possibility of 
non-classical light generatiouj^^"'*^ or single-photon opti- 
cal switching!^ Beyond that, short-distance coupling in 
quantum dot "molecules" has been demonstrated)^"— 
where however the coupling is enforced by the direct 
overlap of the QD wave-functions and/or the electro- 
static Forster dipole-dipole interaction;^!^ rather than 
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by any long-distance mechanism. Altogether, these ad- 
vances suggest that the field has reached the milestone, 
following which the process of long-distance, photon- 
mediated interaction between two or more quantum dots 
should also be addressed. It has been shown that the 
light-matter interaction between a QD and the electro- 
magnetic modes of a non-structured photonic environ- 
ment is very weak^^— thus photonic structures are 
needed in order to tailor the density of optical modes 
and thus enhance radiative coupling between spatially 
separated quantum dots. Indeed, short-range radiative 
coupling has already been achieved in several experi- 
ments involving small optical cavities, where strong cou- 
pling of two quantum dots to the same cavity mode was 
detected,—"— and, most recently, its coherent nature was 
demonstrated)^ A photonic structure has also brought 
the experimental demonstration of long-distance transfer 
of photons emitted by an embedded QD.™ 

On the theoretical side, specific aspects of structures 
with one or more quantum dots in a photonic envi- 
ronment have been studied. These include the strong 
coupling regime and emission spectrum of one^*"^^ or 
more^"— dots in a microcavity, as well as the possibility 
of performing cavity-mediated qubit operations through 
coherent excitation exchange in such a systemf22r2S 
In addition, the spontaneous emission enhancement of 
one dot coupled to a single waveguide mode has been 
estimated.—"— There are, however, only a few studies of 
the dot-dot interaction at a mesoscopic (i.e. more than 
one wavelength) inter-dot distance - which is a main fo- 
cus of this work. Most notably, the possibility to generate 
entangled states between distant QDs in a coupled-cavity 
system was recently demonstrated,^*^ as well as the non- 
trivial decay dynamics^.^ of two distant dots in a photonic 
crystallitci^ However, a general formalism accounting for 
an arbitrary number of quantum dots coupled to arbitrar- 
ily many photonic modes is still lacking, and in particu- 
lar, the distance dependence of the radiative interaction, 
the influence of fabrication disorder, and the competition 
between excitation transfer at-a-distance and radiation 
losses still remain open questions. To address those, a 
microscopic description of light-matter coupling with a 
realistic description of the photonic modes is needed. 

In this paper, we lay down the semi-classical linear re- 
sponse theory for a system of N distinct, spatially local- 
ized excitonic transitions in QDs, coupled to M photonic 
modes of an arbitrary photonic structure. In particular, 
we frame the underlying Maxwell equations into an eigen- 
value problem, describing the polariton modes of the sys- 
tem in analogy with the polariton formalisms for a bulk 
semiconductor—, for quantum wells^, and for QDs in an 
unstructured photonic environment.- 51 Por the compu- 
tation of the photonic modes, the Bloch-mode expansion 
method is employed;^ although any other method which 
provides reliable field profiles (e.g. FEM, FDTD) can also 
be used. Even though the modeling of radiative effects 
in presence of fabrication disorder lies beyond the scope 
of the present work, the Bloch-mode expansion is partic- 



ularly well-suited for treating large, disordered photonic 
structures^i^ and was thus an obvious choice in view 
of a future extension to disordered PHCs. We apply the 
formalism to Ln cavities and W 1 waveguides based on 
a PHC slab. We show how known single-dot radiative 
properties - such as the vacuum Rabi splitting in a mi- 
crocavity and the Purcell enhancement and /3-factor in 
a waveguide - are well reproduced. The main focus of 
the work however is the quantitative characterization of 
radiative coupling between two dots in those same struc- 
tures. To this purpose, we characterize the spectra of the 
polariton eigenmodes, the time-evolution of a starting ex- 
citation in one of the dots, and the distance-dependence 
of the radiative excitation transfer in a spatially extended 
structure. Our simulations provide a comprehensive pic- 
ture of the effective dot-dot radiative coupling, and show 
that, with realistic PHC and QD parameters, a sizable 
interaction can be expected at mesoscopic distance. 

The work is organized as follows. In Section II we 
derive the main theoretical formalism, while in Section 
III we thoroughly discuss the values of the parameters 
entering the model, for realistic InAs/GaAs-based semi- 
conductor nanostructures. Section IV contains the main 
results obtained from the application of the model to the 
study of one and two QDs embedded in Ln cavities and 
W 1 waveguides. In Section V we present our conclusions 
and an outlook of future work. 



II. THEORETICAL FORMALISM 

Starting from Maxwell's equations with the assump- 
tions of a non-magnetic medium and no free charges, the 
electric field in the frequency domain obeys the equation 
(written in Gaussian units) 

^2 

V X V X E(r, uj) (e(r)E(r, uj) + 47rP(r, w)) . (1) 

In particular, here the spatial dependence of the dielec- 
tric constant, e(r), completely characterizes the under- 
lying photonic structure, while the optical response of 
the quantum dots is included in the polarization vector 
through a non-local susceptibility tensor^, such that 

F{r,Lu)= J dr'x(r,r',c^)E(r',w). (2) 

In what follows, we will consider the specific case of 
excitons originating from the heavy-hole band of a semi- 
conductor with cubic symmetry (e.g. InAs), for which 
only the x— and y— components of the polarization cou- 
ple to the electromagnetic field according to the following 
susceptibility tensor— i^i^S 
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The formahsm can be easily generalized to different forms 
of the susceptibility tensor. Here, a runs over all QDs, 
/i^^, is the squared dipole matrix element of the inter- 
band optical transition, (r) — ^^(re = r,r/i = r), and 
^'Q(re,r^) is the excitonic wave- function, normalized as 



drey dr,,|^'c„(re,r,,)|^ 



= 1 



(4) 



We denote the frequencies of the bare excitons by a su- 
perscript a, in order to distinguish them from the fre- 
quencies of the photonic resonances, which we will later 
on index with subscripts, e.g. as uJm- Notice also that 
here all frequencies are assumed to be complex quan- 



tities, e.g. u}'^ 



3fi(a;") 



-, where 7" represents the 



overall decay rate of the exciton state, including any pos- 
sible non-radiative mechanism and the rate of radiative 
decay into photon modes that are not included among 
the M modes treated exactly. 

In order to turn the Maxwell equation into a self- 
adjoint form, we introduce the quantities^'' Q(r, cj) ~ 
\/£(r)E(r, oj). Eq. JT]) then becomes 



The sum in Eq. ([8]) runs in principle over the infinite 
set of eigenmodes. In most situations of interest, how- 
ever, this sum is dominated by the resonant modes of the 
photonic crystal that are closest to the frequency range 
characterizing the excitonic transitions. In addition, in 
all structures of interest (e.g. a PHC^i^^, pillar cavity^ 
or a microdisc*^), the dots are typically embedded within 
the dielectric medium, i.e. their wave-functions are non- 
negligible only in a region where e(r) = £00, the permit- 
tivity of the semiconductor. Thus, as the r-dependence 
of all quantities will eventually enter through overlap in- 
tegrals w ith t he QD wave-functions, we can safely sub- 
stitute \/e(r) = ^/e{r") = ^/e^. Finally, in typical sit- 
uations, all QD transition frequencies lie within a small 
range originating from the inhomogeneous distribution 
of QD sizes. A very good approximation consists then 
in replacing the ui on the r.h.s. of (O, as well as the 
(w,„-|-w)/2 obtained by factoring the denominator in 
with an average exciton transition frequency wq • In order 
to compute the complex frequency poles, corresponding 
to the resonances of the coupled system, we consider the 
homogeneous problem associated with Eq. ([T])- Then, by 
defining 



TQ(r,L.)-— Q(r,L^) = 



J dr'x(r,r',w; 



(5) 



which is an inhomogeneous differential equation defined 
for the self-adjoint differential operator 



we obtain 



Q"(w) = J dr^'„(r)Q(r,w), 



(10) 



T 



1 



=V X V X 



1 



(6) 



The susceptibility tensor as given in Eq. |3] decouples the 
z-polarized fields. We then define the two-dimensional 
field Q = {Qx,Qy)- We can solve the problem using a 
Green's function approach^, in which the formal solution 
to Eq. dl]) is 



Q(r,c^) 



' y^y^ Q„,(r)0Q°r 



Q"(w). 



By integrating Eq. [7] with J dr'I'^(r) and defining addi- 
tionally Q"(w) = Q"(a;)/(w" — w), we finally obtain a 
set of equations (labeled by P) for the complex frequency 
poles 



Q(r,w) = Qo(r,t^) 
47r 



(7) 



dr' / dr"G(r,r-,^) ^^'';f;"^^ Q(r",a.). 

The Green's tensor can be expanded onto the basis 
of field eigenmodes using the resolvent representation, 
following Fredholm's theory^ 



Q— 1 m— 1 



We now define the quantities 



Q"(c.). 
(11) 



G(r, r',tj) 



Q^(r) ® Qr„(rO 



(8) 



where the Qm-s are the orthonormal eigenfunctions of T 
corresponding to eigenvalues us^jc?, and ® is an outer 
product defined as 



A(g) B 



■^x -^x 

AyBx Ay By 



(9) 



/o 2 \ 1/2 

\ £00 h 



Q'L. (12) 



which should be interpreted as the coupling strengths 
between the m-th mode of the PHC and the a-th QD. 
To this end, we notice that the 2N equations in pT|) can 
be solved only for those values of uj for which the N x N 
matrix 
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Ai = 



/ 1 v^M 



\ 



1 V^J« 



EM 
m—1 



UJ^ —UJ 
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m—1 


UJ 


m — 



(13) 



A^m—l 



UJ. 



N 



is singular. This is a nonlinear equation, but we notice 
that it can be transformed into a more familiar form, 
since it is mathematically equivalent to finding the eigen- 
values of the matrix 



for the case of two quantum dots embedded in several of 
the most widely studied photonic crystal structures: the 
L3 and Ln cavities^^, and the Wl waveguide. 







A, = 



\9 



1* 



5aL 



9lx 



9l,y 



N* 
yM,y 



9m,x\ 

9li, 



9M,y 





(14) 



More precisely, solving det(Ai) = is equivalent to solv- 
^I{2NxM)x{2NxM)) = 0, whenever uj ^ 



ing det(A2 

uimyrn = 1...M. The proof can be easily obtained 
by, on one hand, multiplying the equation for Ai by 
Y[m=ii^rn — uj), and on the other, using in the eigenvalue 
problem for A2 the following identity for the determinant 
of a block-matrix: 



det 



A B 
C D 



det{D) det{A~ BD~^C) . 



(15) 



The poles w = uim will generally exist as solutions 
only when a photonic mode Q™ is fully decoupled from 
the system, i.e. when = OVa, in which case this 
mode can safely be excluded from the very beginning. 
The 2N + M complex eigenvalues of A2 then define the 
frequencies (real part) and the loss rates (— 2x imaginary 
part) of the polariton modes of the system, while the 
eigenvectors 



— i-^x^ Kit ■ ■ ^x ^ J -^1 • ■ • : Aj\/) 



(16) 



define the corresponding Hopfield coefficients^^. Notice 
in addition that the matrix of Eq. ([Ti| corresponds to 
a Tavis-Cummings Hamiltonian^^ in the weak excitation 
regime, when only transitions from the ground state to 
the manifold of states with a single excitation are con- 
sidered. 

The present formalism applies to a large variety of pho- 
tonic structures and to an arbitrary spatial distribution of 
QDs. In this sense, it generalizes the results that were ob- 
tained for specific configurations^i^"— . As an illustrat- 
ing application, in section HVl we present results obtained 



III. MODEL PARAMETERS 

In order to quantify the susceptibility ([3|), we need an 
appropriate model of the exciton wave- function evaluated 
at equal electron and hole positions, ^'^(r) = 'i'ai^e = 
r,r/i = r). This function is not properly normalized as a 
function of r (the correct normalization is over R'^ x R'^ as 
given in Eq. (U)). In fact, similarly to the quantum well 
cas o^^'^^ , the oscillator strength of the exciton transition 
in the QD depends on the dimensionless quantity 



dr*„(r) 



(17) 



The particular shape of the wave-function enters 
through the overlap integrals with the electric field, as 
given in Eq. ([TU]). As long as the size of the QDs 
is much smaller than the characteristic wavelength, the 
electric field varies very weakly in the region where ^'q, 
is non-negligible, and thus the point dipole assumption, 
\I'a(r) = C6{r — Fc), is a very good approximation. In 
what follows we will mostly use parameters typical of self- 
organized InGaAs QDs^S., whose size lies in the 10 — 20nm 
range, with a typical exciton recombination energy of 
1.3eV (A « 950nm). For these values, we checked that 
assuming a Gaussian shape for ^'^(r) introduces little 
change with respect to the Dirac-delta assumption. No- 
tice, however, that the strong dependence2ii2^ of the QD 
oscillator strength with its size is still present, carried by 
the normalization constant C. One way to estimate this 
constant is through a microscopic model of ^'^(r)--'^.. 
Here, instead, we take a more pragmatic approach, and 
compute C based on the measured radiative decay rate 
of QDs. Following Ref. [El], this is given by twice the 
imaginary part of the quantity 



'dfc|vl/.^p ^(^^° (18) 



where ^'ck is the Fourier transform of 4'a(r). With the 
assumption vE'a(r) — CS{r — Va), the decay rate is thus 



4- P 

3 hSoo 



(19) 
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where ko = {uJo/c)^/e^, and we defined the dipole mo- 
ment d of the dot (also labeled or p*^) as 



(20) 



Eq. ()19p coincides with the expression that is com- 
monly adopted^i^saoi. For typical QDsi^ii^JH^ with 
radiative lifetime of Ins and cxciton transition energy 
fiw" « 1.3eV, we obtain a squared dipole moment cP w 
O.SleV X nm^. 

The last requirement of the problem is the knowledge 
of the modes of the PHC structure, i.e. the set of or- 
thonormal functions {Q,„(r)} and their corresponding 
eigenfrequencies ujm- Here, PHC modes are computed 
using the Bloch-mode expansion method^^, which con- 
sists in expanding the modes on the basis of the Bloch 
modes of a regular waveguide. These latter were in turn 
computed using an expansion over the guided modes of 
a uniform dielectric slabJ^ This approach turns out to 
be particularly well suited for elongated PHC cavities 
as considered in the present work. The computation 
was carried out over a finite supercell S in the plane 
of the crystal, and infinite space along the orthogonal, 
z— direction. The orthogonality relation is then given by 



/ d^p r dz Q„(p,z)Q:(p,z) 

J S <J —oo 



(21) 



All the photonic crystals we consider are based on a 
triangular lattice of circular holes etched in a dielectric 
slab suspended in air. The specific parameters we chose 
are relevant to GaAs structures j^i^iiS^ namely: lattice 
constant a = 260nm, hole radius 65nm, and slab thick- 
ness 120nm, with a permittivity ^/e^ = 3.41. In this 
work we consider only ideal PHC structures, in the ab- 
sence of any disorder that would arise from the fabrica- 
tion process. Disorder in PHCs has two important ef- 
fects. First, it determines extrinsic radiation loss rates of 
otherwise fully guided modes in waveguides, and strongly 
suppresses the quality factors of high-quality PHC cav- 
ities. This effect is here taken into account through the 
inclusion of a constant phenomenological loss rate for the 
modes under study, related to their quality factor Q by 
7 = i^/Q- For the L3 cavity of section HV Al we set 
Q = 10000 or 30000. For the longer Ln cavities (sec- 
tion [IVB]) and the Wl waveguide (section llV Cp . we set 
Q = 50000 for all modes. The second way disorder af- 
fects the results is by modifying the spatial profiles of the 
electric field modes, especially in the case of waveguides. 
This effect lies beyond the scope of the present work - 
although a brief discussion is given in section |V] - and 
will be the object of a future work. 




A. Application to an L3 cavity 

The system of one quantum dot coupled to an L3 cav- 
ity has been widely studie d'^^i"^^'^"^ and is thus a good 
starting point for testing the present formalism. 
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FIG. 1. (Color online) Imaginary part of the electric field of 
the fundamental mode of an L3 cavity, (a); 'S{Qx{r)) and 
(b): 9((5y(r)). In the one-QD simulation, the dot was placed 
in the central maximum of the y-field (dot position marked 
by a white cross). For the two-QD simulations, the dots 
were placed in the corresponding secondary maxima (posi- 
tions marked by white stars). 

The cavity is a modified L3 cavit y^^'^"^ , where the two 
holes on each side of the cavity are shifted outwards by 
0.15a, and their radii are decreased by 80%. This de- 
sign improves the quality factor by more than one decade 
compared to that of a standard L3 cavity, while chang- 
ing the field profile only marginally. We include in the 
computation only the fundamental cavity mode, shown 
in Fig. [H We further assume the QD to lie on the in- 
plane symmetry axis of the cavity, where Qx = 0. The 
diagonalization of the matrix (jl4p is then equivalent to 
the well-known expression 



detr^T'^ 9^ UO. 
V 9c ^c-i^J 

The coupling constant gc, through Eq. is 



(22) 



9c 



dQyiVa) ; 



(23) 



which matches previous theoretical results^ when the 
dot is sitting in the center Tq of the cavity and the mode 

|(5j,(ro)p. As expected from 
- 7yp/16, vacuum- field Rabi 
splitting appears between two polariton modes. The en- 
ergy splitting at zero dot-cavity detuning is given by 2?if2, 
where the Rabi frequency n is 



volume is defined as y = 
Eq. for |5cP > |7. 



IV. APPLICATIONS 

In this section, we apply the formalism to the proto- 
typical cases of one or two QDs embedded in elongated 
Ln cavities or in a Wl waveguide. 



9c\ 



(7c " 7.;)^ 
16 



(24) 



Using the PHC and QD parameters we already intro- 
duced, the coupling constant was computed to be h\gc\ = 
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147/^eV, which compares perfectly with the most recently 
reported result for that system^. 

After showing the way the standard one-dot cavity- 
QED results are reproduced with our formalism, we now 
proceed to the situation of two dots coupled to the same 
cavity mode (see Fig. [T]), and so radiatively coupled to 
each other. 

We assume a symmetric spatial configuration of the 
two dots with respect to the cavity center (see Fig. [T|), 
resulting in equal coupling constants hlgd = 125/ieV 
for the two dots. Since usually 7c ^ jy, i.e. the 
losses through the cavity mode are significantly larger 
than the QD losses through other channels (both non- 
radiative and radiative through modes other than the 
cavity mode), we set here and in all following sections 
7^'^ = 0. Given the phenomenological way these rates 
enter the formalism, calculations can easily be general- 
ized to include finite QD loss rates. Let us first consider 
the case of zero dot-dot detuning S = ojy — ojy. The rel- 
evant exciton states are in this case the symmetric and 
antisymmetric linear combinations of the two QD states, 
whose coupling to the cavity mode depends on the sym- 
metry of the elect ric field profile. As discussed exten- 
sively in Ref. |l04l | , the L3 cavity symmetry is described 
by the D2h point group, and its fundamental mode be- 
longs to the B2u irreducible representation, which is even 
with respect to the ayz symmetry operation (as can also 
be seen from Fig. [1]). Hence, the antisymmetric QD 
state remains dark, while the symmetric one behaves as 
a single exciton with a coupling constant y/2gc- 

In Fig. [2] (a) we plot the eigenfrequencies of the system 
as a function of the detuning between the exciton reso- 
nance frequency ojy (same for both dots) and the cavity 
resonance frequency ujc, as computed for a cavity quality 
factor of Q = 10000. We observe vacuum Rabi split- 
ting between an upper and a lower polariton in exactly 
the same way we would for a single dot coupled to the 
cavity, but in addition we see a dark mode which is a 
trivial solution, uj — Uy. The splitting between the lower 
and the upper polaritons at zero dot-cavity detuning is 
2Mlc — 347/ieV, which for Q = 10000 corresponds ex- 
actly to an effective coupling constant of -\/2 x 125/^eV. 
The system is further characterized in panels (b) and (c) , 
where we plot the Hopfield coefficients for each of the 
three eigenmodes (correspondingly color-coded). This 
clear collective behavior has been observed experimen- 
tally in a QD-cavity systent^"— , while the more general 
dependence of the effective coupling constant with the 
number of coupled two- level systems N - given by -s/^lffcl 
- has also been observed in a circuit-QED system It 
is very important to remark that this dependence has 
nothing to do with the ^/N\gc\ energy splitting of differ- 
ent rungs in a Jaynes-Cummings model, where N would 
be the number of photons in the system: on the contrary, 
as discussed before, here we restrict to the linear response 
only, which holds in the limit of vacuum electromagnetic 
field. The effect in our case is simply due to the collective 
behavior of the TV resonant quantum dots. 




-1000 -500 500 1000 



FIG. 2. (Color online) (a): Eigenfrequencies (solid lines) 
and radiative rates (dashed lines) for two QDs with no dot- 
dot detuning, strongly coupled to an L3 cavity mode with 
Q — 10000. With a dashed-dotted line, the bare cavity res- 
onance is also indicated. The Hopfield coefficients for each 
solution, correspondingly color-coded, are presented in panels 
(b): equal (in absolute value) QD coefficients and (c): cavity 
coefficient. 



The major experimental challenges to the radiative 
coupling of two spatially separated quantum dots is 
achieving both spatial control (to ensure strong overlap 
between each of the dots and the cavity mode) and spec- 
tral control (to ensure as small dot-dot and dot-cavity de- 
tuning as possible). Typically, QDs are characterized by 
an inhomogeneous distribution of exciton energies with 
a width of several meV. Then, two QDs are very likely 
to be detuned. In Fig. [31 we study the same system, 
but assuming a detuning S — 300/ieV. Close to reso- 
nance, all of the eigenmodes acquire a finite component 
from the cavity mode. Additionally, they have both a 
significant |A^| coefficient (panel (b)), and a significant 
|A^| coefficient (panel (c)), implying that there is a siz- 
able radiative coupling present. The radiative coupling 
is expected to vanish as the cavity-dot detunings become 
much larger than the coupling constant, and an expres- 
sion for an effective coupling strength in this limit was 
derived in Ref. [zl and il06| . 

We now address the question of how the excitation 
transfer process depends on time. This aspect is of partic- 
ular importance to assess the usefulness of the radiative 
excitation transfer as a coupling mechanism between dif- 
ferent qubits in a semiconductor-based quantum gate ar- 
chitecture. In the present case, when polaritonic features 
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FIG. 3. (Color online) (a): Eigenfrequencies (solid lines) and 
radiative rates (dashed lines) for two QDs with a dot-dot de- 
tuning of 300/ieV, strongly coupled to an L3 cavity mode 
with Q = 10000. With dashed-dotted lines, the bare excitons 
and the bare cavity resonances are also shown. The Hopfield 
coefficients for each solution, correspondingly color-coded, are 
presented in panels (b): first exciton coefficient, (c): second 
exciton coefficient, and (d): cavity coefficient. 



are spectrally resolved, one correspondingly expects the 
excitation to oscillate between the different basis states, 
including the photon state. To illustrate this aspect, we 
compute the time-dependent amplitudes of the various 
basis states, assuming that one QD is excited at t = 0. 
From these amplitudes, we extract time-dependent prob- 
abilities of finding the excitation in each of the basis 
modes, expressed in vector form as 



m = 



A, 



(25) 



where A2 is the matrix of Eq. ([H)) . These probabilities 
are properly normalized if one accounts also for the prob- 
ability Poutit) of the excitation to have radiated out of 
the system, i.e. ^Pi{t) = 1 — Pout{t)- In Fig- IHwe plot 
these time-resolved probabilities for a starting excitation 
in one of the QDs, i.e. Ai„ = (1,0,0). We study four 
different cases: either zero dot-dot and dot-cavity detun- 
ing, or h6 — 300/ieV (with the cavity frequency tuned 
at the average of the two exciton frequencies), and cav- 
ity Q-factor equal to either 10000 or 30000. In panels 
(a) and (b), where 6 = 0, the probabilities never decay 
to zero due to the presence of a dark state and the fact 
that no non-radiative decay mechanism was included. In 
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FIG. 4. (color online) Time evolution of the probability of an 
excitation in one dot to be transferred to the second dot or 
to the cavity, (a): 5 = 0, Q = 10000; (b): 5 = 0, Q = 30000; 
(c) hS = 300AteV, Q = 10000; (d) M = SOO/xeV, Q = 30000. 



panels (c) and (d) a dark state no longer exists, and a 
clear decay of the excitation with a characteristic life- 
time depending on the Q-factor is visible. All plots show 
that the excitation oscillates between the three possible 
states, on a time scale defined through the radiative cou- 
pling strength. In particular, the probability of finding 
the system in an excited state of the second QD remains 
sizable over several oscillation periods, showing that a 
significant dot-dot interaction can be achieved with ex- 
perimentally feasible parameters. These results generally 
agree with specific setups of radiatively coupled QDs in 
photonic crystals, that have been recently studied in the 
literature^SiS^ 



B. Application to Ln cavities 

Recently, using Ln cavities with n > 3 to achieve light- 
matter coupling has spurred interest^iSl^iSS, as these cav- 
ities generally have a larger quality factor than the L3 
ones - though at the expense of a larger mode volume 
and thus a smaller dot-cavity coupling strength. In this 
section, we take the Lll cavity as a first case study, since 
coupling of quantum dots to its photonic modes has al- 
ready been experimentally demonstratedi^liiSS. 

Fig. [5] shows the first four modes, Mi_4, 
of the Lll cavity, with resonant energies 
1.3065eV, 1.3125eV,1.3269eV, and 1.3565eV, re- 
spectively. The two dots were placed in the center of an 
elementary cell on each side of the center of the defect 
(i.e. at a distance a from the center of the cavity and 
so 2a from each other), where the coupling constants 
for each of them are \hgi\ = 94/zeV, \fig2\ = 55/^eV, 
|%3| = 65/zeV, and \hg4\ = 89^eV. Since the smallest 
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FIG. 5. (Color online) Q(Qy{r)) for the four lowest-energy 
modes of the Lll cavity; (a): Fundamental mode Mi, at 
hcoi = 1.3065eV, (b): Ma, fe^2 = 1.3125eV, (c): M3, riios = 
1.3269eV, (d): M4, &J4 = 1.3565eV. The positions of the 
quantum dots are marked with white stars. 



energy difference between the cavity resonances is 
between uji and lu2, and is ~ 6meV, i.e. much larger 
than all the coupling strengths, it is reasonable to expect 
that the dots will never couple significantly to more than 
one mode. 
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FIG. 6. (Color online) (a): Eigenfrequencies for 2 QDs with 
no dot-dot detuning with exciton resonant frequency in the 
range of the two lowest cavity modes (indicated by dashed- 
dotted lines) of an Lll cavity, with quality factor of each 
mode Q = 50000. The Hopfield coefficients for each solu- 
tion, correspondingly color-coded, are presented in panels (b): 
equal (in absolute value) QD coefficients, (c): coefficient of 
Ml, and (d): coefficient of M2. 



This is indeed what we observe in Fig. |6l where we 
plot the eigenfrequencies and Hopfield coefficients of the 
system with zero dot-dot detuning, as a function of the 
exciton transition frequency ujy — ujy. The Q value for 
all the modes, and for all modes of all Ln cavities consid- 
ered below, was set to 50000 (besides from influencing the 
transition between weak and strong coupling, this value 
has no qualitative effect on the results). In panel (a), 
four distinct polaritons, labeled pi-4, are visible, which 
can again be understood through the symmetry of the 
two cavity modes. The mode Mi (Fig. [3 (a)) is sym- 
metric w.r.t. the ayz reflection, while M2 (Fig. [SJ (b)) 
is antisymmetric. Thus, the symmetric combination of 
the QD excitons couples to Mi but not to M2, and for 
LOy « wi results in two anti-crossing solutions {pi and P2), 
while for « W2 - it gives rise to (p2)- Vice versa, the 
antisymmetric combination of the two QD states is dark 
for LOy uji (ps), and gives two anti-crossing branches at 
u!y w UJ2 {p3 andp4). Just as in the L3 cavity case, when 
some dot-dot detuning 6 is introduced, there are no more 
dark modes, but instead a third polariton eigenmode ap- 
pears in between the lower and the upper solutions. This 
case is presented in Fig. [7] for M = 300/ieV. In both 
cases, this behavior is reproduced also by the variation 
of the Hopfield coefficients for all modes, as displayed in 
the panels below panel (a). 

The situation should change significantly when increas- 
ing the length n of the photonic defect. Then, we expect 
the energy spacing between the resonant frequencies of 
the Ln cavity to decrease and eventually become compa- 
rable to the typical coupling strength. In this situation, 
the radiative transfer process is no longer mediated by 
an isolated cavity mode, and a smooth transition to a 
multi-mode coupling regime is expected. In order to de- 
termine at which cavity length this crossover occurs, one 
should also consider the fact that the coupling of a dot 
to each individual mode decreases with the increase of 
the mode volume. As a result, the crossover length is in- 
creased with respect to what would be given by a simple 
assumption of constant coupling strength per mode. 

In Fig. [SI we plot the minimum mode-separation uj2 — 
LOi VS. the length n of the cavity, and in addition show 
the coupling strengths for m = 1 ... 10. For all n, 
the dots were placed as in Fig. [5]- at a distance a on each 
side of the center of the cavity. The fact that half of the 
coupling constants decay much faster as a function of n, 
is again explained by the particular symmetry of the field 
profiles. It turns out that for every n, the modes alternate 
between symmetric and antisymmetric w.r.t. ayz, as can 
be seen in Fig. [5] for the Lll case. In the limit of large 
n, the antisymmetric modes have a small amplitude at 
the QD positions close to the node, resulting in a small 
radiative coupling strength. 

The crossover from a single-mode to a many-mode 
regime occurs around n — 141, as clearly visible in Fig. 
m In panels (b)-(g), we show the corresponding Hop- 
field coefficients for three different values of n, given by 
n = 71, n = 141, and n = 211, also indicated by dashed 
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FIG. 7. (Color online) (a): Eigenfrequencies for 2 QDs with 
dot-dot detuning h& — SOOpeV in an Lll cavity, Q = 50000. 
The corresponding Hopfield coefficients are also presented, 
(b): first dot, (c): second dot, (d): of M\, and (d): of M2. 



lines in panel (a), and for two different polariton modes. 
Consequently, for n = 71, the Hopfield coefficients of pi 
and p2, shown in panels (b) and (e) respectively, are still 
largely dominated by one cavity and one dot component. 
On the other hand, for n = 211 - panels (d) and (g) - the 
value of several photonic fractions Am is non-negligible. 

In Fig. HI we plot the polariton energies as a function 
of QD-exciton energy in the case of the LUl cavity, for 
(5 = 0. As mentioned already, the photon modes alter- 
nate between symmetric and anti-symmetric w.r.t. the 
^yz operator, hence coupling to either the symmetric or 
anti-symmetric linear combination of the QD states is 
present. In the figure, the polaritons due to a combina- 
tion of symmetric states are denoted by blue lines, while 
the anti-symmetric combinations are represented by red 
lines. In the symmetric case, the exciton-photon coupling 
strength is always large, and anti-crossing occurs at ev- 
ery mode. In the anti-symmetric case, the results show 
a transition from weak coupling (close to the lowest i^ra) 
to strong coupling (anti-crossing is visible in the higher-w 
inset), due to the fact that the coupling strengths there 
become larger than uJml^Q,- It is clear both from Figs. 
|S] and ini that for n — )■ 00, the dots couple to a structured 
continuum of photon modes, reproducing the physics of 
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FIG. 8. (Color online) (a); Black lines - coupling constants 
between one QD and the ten lowest modes of an Ln cavity, 
vs. n; green line - energy separation between the lowest two 
cavity modes, (b)-(d): Hopfield coefficients of pi (see Fig. |6]) 
as a function of the bare exciton frequency uj^ with no dot-dot 
detuning, for n = 71, n — 141, n — 211 (the values marked 
by dashed vertical lines in (a)). The red line shows the dot 
coefficients, while the blue lines belong to the many cavity 
modes, (e)-(g): Same as (b)-(d) but for the polariton p2. 



a Wl waveguide. This situation is studied in the next 
section. 



C. Application to a Wl waveguide 

The results obtained for a Ln cavity indicate that ra- 
diative coupling is still sizable in very long structures and 
might be effective even at very long distance between the 
two QDs. Here, we investigate this possibility in more 
detail, by considering QDs embedded in a Wl photonic 
crystal waveguide. 

Coupling of a single dot to a Wl waveguide (or a sim- 
ilar structure) with the purpose of spontaneous emis- 
sion enhancement (and the potential application as a 
single-photon source) has already been widely discussed 
theoreticallyiS^iii^, and achieved experimentallj^"— . 
The fact that it is already possible to couple efficiently a 
dot to the guided modes of the waveguide is promising 
in view of achieving radiative coupling between two dots 
that could - due to the spatial extension of the struc- 
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FIG. 9. (Color online) Eigenfrequencies for two QDs with 
no dot-dot detuning in an LI4I cavity (Q = 50000 for each 
mode) vs. the resonant frequency of the excitons. The insets 
show close-ups over two selected regions. 



tures and the modes they support - extend to inter-dot 
distances for which targeting each dot individually by a 
laser pulse is possible. 

We begin our study by looking at the modes of the Wl 
waveguide. The band structure is presented in Fig. [TU] 
(a), where two guided bands in the band-gap of the reg- 
ular crystal are visible. Strongest dot-PHC coupling is 
typically achieved for the smallest group velocity (largest 
local density of states of the photonic modes), and so the 
spectral range we concentrate on is around the band-edge 
of the main guided band (blue line) , where the group ve- 
locity of the ideal photonic structure vanishes. The sec- 
ond guided band is spatially odd with respect to a axz 
reflection, ^''■^ and would not couple to the exciton state 
of a QD located at the center of the waveguide. In the 
simulations below, we compute the Wl modes for 2048 
fc-points in the interval (— 7r/a, 7r/a], which is equivalent 
to simulating a waveguide of length 2048 elementary cells 
with periodic boundary conditions (PBC). In panels (b)- 
(e) of Fig. [TUl we show the electric field profiles of four 
modes lying close in energy to the band edge of the main 
guided band. As is the case with all structures we consid- 
ered so far, this band has vanishing component on the 
symmetry axis of the waveguide, allowing us again to in- 
clude the t/-polarized fields only. Furthermore, modes at 
±fc are degenerate - one propagating and one counter- 
propagating - with real-space profiles proportional to 
exp{ikx) and exp(— jfca:), respectively. As basis states, 
we take the symmetric and the anti-symmetric combi- 
nation of the degenerate guided modes, representing the 
fields by their "standing wave" profiles: one with a max- 
imum and one with a zero amplitude in the center of the 
guide (compare panels (b)-(c)). Without loss of general- 
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FIG. 10. (Color online) (a): Band structure of the Wl waveg- 
uide; the dashed line shows the light cone. The QD resonant 
energies are close to the band-edge energy of the main guided 
band (blue). The field profiles of four guided modes in that 
spectral region are shown, over a small stretch of the waveg- 
uide: (b), (c): the two degenerate modes at hojk ~ 1.30308eV; 
(d): the symmetric mode at hwk ~ 1.30224eV (e): the sym- 
metric mode at hiJk = 1.30218eV (at the band edge). In all 
computations, one QD was placed in the center of the waveg- 
uide (white cross), while the second one - in the center of one 
of the successive elementary cells (white stars). 



ity (due to the PBC), we place one dot at that position, 
so that it couples to one of the modes only. The second 
dot is then placed in the center of a successive elementary 
cell, and will, in general, couple to every mode in the ba- 
sis thus constructed, so even for zero dot-dot detuning, 
no fully dark state is present. 

In Fig. [11] (a), we show the polariton structure in 
the spectral range close to the band edge of two dots 
with dot-dot detuning S ~ Q, with the second dot placed 
at the closest possible distance, a, from the first one. 
The quality factor of each of the photonic modes was 
again set to Q = 50000 for all modes. While a strong 
dependence of the Wl loss rates on the group velocity 
close to the band edge has been shown in transmission 
measurementsiii, this dependence is heavily influenced 
by back-scattering due to disorder. In our case, we model 
stationary modes rather than transport, and the only rel- 
evant radiative loss is the one out of the plane of the 
PHC slab. Then, the assumption of approximately con- 
stant Q-s is realistic, as seen from microscopic model- 
ing of extrinsic disorder-induced losses^^^i^. Polariton 
modes originating from antisymmetric photon modes are 
essentially uncoupled and are not displayed (although, 
they were still included in the computation). The cou- 
pling constants of each of the dots to each of the symmet- 
ric modes varies very little, and is fi\gli^\ ~ 7/ieV. The 
uj = ujy solution (straight diagonal in panel (a)) is due 
to the anti-symmetric QD combination, which is almost 
dark. The strongest anti-crossing behavior is exhibited 
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FIG. 11. (Color online) (a): Eigenfrequencies for two QDs 
in the waveguide, with no dot-dot detuning. The Hopfield 
coefficients for the red line of (a) are shown in (b), where 
the green line shows the QD coefficients |Ay| = |Ay[, and the 
black lines, all the waveguide coefficients |Am|. Same in (c), 
but for the blue line of (a). The loss rate of each mode is 
7w = ajJ/50000 « 26^ieV. 



by the polariton lying below the band edge (blue line), 
whose Hopfield coefficients are given in panel (c). The 
remaining polariton modes display similar behavior, so 
the Hopfield coefficients of just one of them (the red line 
of (a)) are given in panel (b). For completeness, the 
same plots but for M = 100/ieV are given in Fig. [121 
In this case, no dark modes are present and Hopfield 
coefficient corresponding to the two QDs are generally 
different from each other, as seen in panels (b) and (c). 
In both Figs. [TT] and [12] we observe that anti-crossings 
are still present - though characterized by a very small 
energy splitting - where the exciton becomes resonant 
with the various guided modes. This situation can be 
understood as the precursor to the structured continu- 
ous spectrum of modes that would arise in the limit of 
infinite waveguide length, analogous, for example, to the 
polariton modes arising from the interaction between an 
exciton in a two-dimensional quantum well and the three- 
dimensional continuum of electromagnetic modes i^^'^^ 

The present formalism provides a detailed quantitative 
account of the effect of the guided electromagnetic field 
on the radiation properties of few QDs. In particular, 
we derive below the Purcell enhancement of the radia- 
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FIG. 12. (Color online) (a): Eigenfrequencies for two QDs 
in the waveguide, with dot-dot detuning h5 = 100/^eV and 
7w as in Fig. 1111 The Hopfield coefficients for the red line 
of (a) are shown in (b), where the two green lines show the 
QD coefficients lAj''^], and the black lines, all the waveguide 
coefficients |Am|. Same in (c), but for the blue line of (a). 



tive rate characterizing a single QD, and the distance- 
dependence of the radiative excitation transfer process 
between two distant QDs. We compute these properties 
both numerically, and analytically. To this purpose, let 
us consider the elements of the matrix Ai introduced by 
Eq. dH: 



A 



afi 



(t^"-c^)<5„^-G"'^(w), 



(26) 



where the coupling matrix elements G"'^ are proportional 
to the Green's function of eq. [H 



M 
in—l 



.9m.% 



27r uj 



2 



■G{v^,vp,u). (27) 



For a structure with no sharp resonances - like the waveg- 
uide - we can take advantage of the exciton-pole ap- 
proximation and substitute lj — ujq in the denominator, 
in which case G'""(wo) is the self- interaction energy of 
each dot, while |G'^^(a;o)| is an effective coupling constant 
for the case of two dots with zero dot-dot detuning, i.e. 
ujy — — ^0- In order to derive an analytical expression 
for the coupling, let us replace the sum with an integral 
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over k = kx, and use the fact that, in accordance with 
Bloch's theorem, when and are in the center of an 
elementary ceh, gki^p) = exp {—\kx)gk{Ya), to write 



dfc 



\9k\ 



27r J _E. uj{k) — cjo 



(28) 



A few simphfications are due. First, we write ujk — 
lk{(jj{k)) and 7w — — 29(cij(fc)), and assume the latter 
is constant, equal to ujq/Q. Furthermore, we assume 
Ifffel^ = |(7|^, i.e. the coupling strength has weak depen- 
dence on k. This feature is due to the small spatial ex- 
tension of the exciton wave function, resulting in a very 
broad distribution in Fourier space with approximately 
constant overlap with all guided modes, and is also con- 
firmed by our numerical results. Finally, by taking fco as 
the positive Bloch momentum for which the guided mode 
is resonant with the exciton frequency wq, and defining 
the corresponding group velocity 



'dk 



(29) 



we get 



a \9\ 



dk- 

^ 



cos (kx) 



ko 



(30) 



2v„ 



This expression holds in the limit where the result- 
ing spectral linewidth is small enough so that the group 
velocity is still well defined. It can now be applied for 

I 



example to obtain the radiative lifetime of a single dot 
embedded in the waveguide as F" = 25 (A""), and so 



tan 



2{k ~ ko)vg 
7w 



k—7r/a 



(31) 



The Purcell factor for the enhancement of the single-dot 
spontaneous emission rate is then given by the ratio be- 
tween Eq. ([3T|) and Eq. (fTQ]). This result takes into 
account the detailed structure of the photonic environ- 
ment resulting from the waveguide. In this respect, it 
generalizes the result obtained by assuming that only one 
Bloch mode at wave vector k = ko determines the radi- 
ation loss process.— 1^ This simplified result is recovered 
by taking the limit 7w — >■ in the integral (P7|. namely 
by assuming that the guided Bloch mode has vanishing 
extrinsic radiation loss rate. The emission rate F' of the 
dot into leaky modes can also be estimated numerically 
by restricting the summation in Eq. (1271) to the modes 
which lie above the light-cone only. Then, the /3-factor in 
the absence of non-radiative decay mechanisms can also 
be computed as 



(32) 



and a further generalization to the case in which non- 
radiative processes are also present follows straightfor- 
wardly. 

Let us now consider the radiative excitation transfer 
process between two QDs. The closed- form expression 
for the cross-coupling term G^^, obtained by carrying 
out the integral ((30|) reads 



^12 ^ a IsT 



cosh I ifco^; ) Ci ( —i— 1- (fc — fco)^; 

fl2 / V ^12 



i sinh I -^^ ikox ) Si ( i (^ ^ ^0)2; 

ri2 J \ri2 



k='- 



k=0 



(33) 



where Ci(z) and Si(z) are respectively the cosine in- 
tegral and the sine integral functions, and we defined 
ri2 = 2ug/7w. The quantity ri2 is simply the decay 
length associated to the propagation of light along the 
resonant guided mode. We expect this decay to char- 
acterize also the distance dependence of the radiation 
transfer process. Indeed, under the ideal assumption of 
vanishing radiation loss rate for the guided mode, in a 
one-dimensional geometry one would expect the radia- 
tive transfer process to be independent of the distance. 
For comparison, as has already been shown, the coupling 
strength decays as in 3D bulk semiconductor--, and 

as R^^J"^ in a 2D planar cavity system^. In Fig. [T51 
we display the absolute value of G^^ computed numer- 
ically through Eq. p7)). for four different values of the 
exciton frequency wq of the two QDs, in a waveguide 



I 

of length 2048a. This quantity is compared to the re- 
sult obtained from the analytical model of Eq. ([33)) and 
to the simpler assumption of an exponential dependence 



IGi^l = IGi 



-2;/ri2 



In panel (d), where wq is taken to 



lie below the edge of the guided band, the group veloc- 
ity cannot be properly defined, and thus the analytical 
model does not apply. 

Apart from this case, it is clear that the distance de- 
pendence of the inter-dot coupling is perfectly captured 
by the simple exponential decay model. The oscillations 
of the numerical curve in panel (a) are due to the finite 
length of the waveguide and reproduce the spatial be- 
havior of the Bloch mode at fc = feg that dominates the 
transfer process. These oscillations cannot obviously be 
reproduced by the analytical model that implicitly as- 
sumes an infinitely extended waveguide. As anticipated, 
the numerical results show that the distance dependence 
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FIG. 13. (Color online) Absolute value of the off-diagonal 
term of the matrix in Eq. (|13p . in the exciton-pole approxima- 
tion, computed numerically for a finite-size waveguide (solid 
Une), analytically through Eq. (|33p (dashed-dotted line), and 
through an exponential decay model with characteristic dis- 
tance ri2 = 2wg/7w (dashed line, 7™ as in Fig. [TT|) . In (a): 
&J0 = 1.30353eV, Ug = 74, (b): huo = 1.30240eV, Ug = 195, 
(c): hojo = 1.30220eV, Ug = 525, (d): hujo = 1.30208eV, i.e. 
100 ^eV below the band edge. 



FIG. 14. (Color online) Time evolution of the probability 
of an excitation in one dot to be transfered to the second 
dot or to the many PHC modes (blue lines). The dot-dot 
detuning is hS = 100/ieV. Horizontally across the panels, the 
inter-dot distance changes from 260nm to 2.6/im to 5.2/im. 
Vertically across the panels, the exciton frequency of the first 
dot changes from (a)-(c): hivl — 1.30224eV (close to the band 
edge), through (d)-(f): huuk ~ 1.30218eV (at the band edge), 
to (g)-(i): fiujk = 1.30208eV (lOO/ieV below the band edge 
energy, which is then resonant with tkjy). 



of the transfer rate is expressed by the decay associated 
to the light propagation, and quantified by the decay 
length ri2. It is interesting to note that even for very 
small group velocities, e.g. Vg < c/500, the interaction 
distance is still of the order of 100a = 26/im, i.e. of meso- 
scopic scale, thus confirming the potential of the Wl for 
very long-distance dot-dot coupling. 

More generally, Eq. (1551) suggests that there is a com- 
promise, enforced by the group velocity, between strength 
and distance dependence of the transfer process. The 
overall strength of the transfer rate depends inversely on 
the group velocity. This expresses the magnitude of the 
local density of states at the QD exciton frequency or, in a 
more suggestive picture, the fact that slow light interacts 
with a QD over a longer time lapse. However, a smaller 
group velocity also implies a shorter characteristic decay 
length ri2, as we are assuming a constant radiation loss 
rate. In a realistic system, including disorder, we further 
expect the group velocity picture to break down at fre- 
quencies close to the band edge, where disorder-induced 
localization of light dominates and the spatial decay as- 
sociated to the localization length becomes shorter than 
ri2. This calls for an analysis including disorder effects, 
that we will consider in a future work. 

We conclude this section by studying the time- 
dependent probability amplitudes of the excitation lying 
in each mode. These quantities are plotted in Fig. [141 
assuming that one QD is excited at i = 0, for three differ- 
ent inter-dot distances and three different values of Uy. 
In all cases, M = lOO^eV was imposed. As discussed 



above, the transfer mechanism is driven by several light 
modes. The plots show that the radiative transfer pro- 
cess still occurs and, in particular, the marked oscillations 
are characterized on average by a period that can be as- 
sociated to an effective transfer rate hfl = 50 — 60/ieV. 
As in short Ln cavities, this rate is quite sizable and 
should be observable in state-of-the-art GaAs-based pho- 
tonic structures 



V. DISCUSSION AND OUTLOOK 

We have developed a general formalism of linear 
radiation-matter coupling in systems of many QDs em- 
bedded in a photonic crystal structure. The formalism 
is an extension of the exciton-polariton formalism well 
known for bulk semiconductors and quantum wells. It 
provides a quantitative account of a variety of radiative 
effects, starting from the basic microscopic parameters 
of the QD-PHC system. As examples of application, the 
main results presented in this work concerned radiative 
effects in the systems of one or two QDs embedded in 
Ln cavities and the Wl waveguide. In the case of one 
QD, we recover the known results for the Purcell en- 
hancement of the radiative rate and the vacuum Rabi 
splitting in the strong coupling regime. In the two-QD 
case, we quantify the strength of the radiative excita- 
tion transfer between spatially separated QDs, which lies 
in the 100 fj,eV range at short distance. The compar- 
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ison of the single-mode coupling strength and the en- 
ergy spacing between modes in Ln cavities of increasing 
length clearly shows that a crossover occurs - around 
n = 150 for GaAs-based systems - between single-mode 
and multi-mode radiative coupling. In the multi-mode 
case, the radiative coupling strength through each pho- 
tonic mode is smaller but the overall effective excitation 
transfer rate still ranges at about 50 /ieV, thus suggest- 
ing that the W 1 is an ideal structure for the realization 
of long-range radiative dot-dot coupling. 

These results suggest that the QD-PHC system could 
be a candidate system to operate as a quantum bus and 
achieve controlled entangling interaction between distant 
qubits. This perspective is corroborated by the two fol- 
lowing remarks. First, semiconductor QDs have recently 
seen a tremendous progress^'^"^^ towards the physical 
implementation of qubits that rely on the electron or 
hole spin as the computational degree of freedom, and 
on the interband optical transition as the main handle 
for single-qubit operations. Second, the optical quan- 
tum bus technology has already been successfully applied 
to achieve controlled two-qubit operations in the system 
of superconducting qubitsji The controlled operation in 
that case has been achieved by moving in and out of 
the anti-crossing region in the polariton spectrum aris- 
ing from radiation-matter coupling. In view of a simi- 
lar development in the semiconductor QD case, at least 
three steps are still needed. First, the ability to fabri- 
cate site-controlled QDs, in order to position them with 
respect to the PHC structure. This is nowadays pos- 
sible thanks to various kinds of growth on a patterned 
substratei^^-iii Second, a clear experimental proof of 
the radiative excitation transfer mechanism at long dis- 
tance, that might only come from ad-hoc technique such 
as, for example, the single-QD two-dimensional four- 
wave-mixing spectroscopy;^ Third, a reliable scheme for 
dynamically controlling the exciton-photon detuning at 
sufhciently high speed. For this latter task, extremely 
promising results are already available on the optical con- 
trol of the resonant frequency of high-Q cavities, partic- 
ularly using carrier-induced optical nonlinearityi^ii^ 



An additional challenge is represented by the task 
of understanding and optimizing disorder effects in the 
light propagation throughout PHC structures. Apart 
from small variations of the field profiles that should 
only marginally influence the magnitude of the cou- 
pling between one QD and a photon mode, the main 
effect of disorder is the localization of light in long 
PHC structures. For the line defects studied in this 
work in particular, localization is known to compete 
with the ability of the waveguide to support slow-light 
propagationjiiSii^ and thus with the enhancement of ra- 
diative effects expected in these structures. More specif- 
ically, when approaching the band edge in a Wl waveg- 
uide, the localization length becomes shorter than the 
decay length related to extrinsic radiation losses. In this 
limit, there is no more light propagation and the concept 
of a group velocity of light no longer holds. As an illustra- 
tive example of this dramatic effect on light propagation, 
the statistical fiuctuations of the transmission through a 
waveguide in the light-localization regime increase and 
the transmission coefficient of a finite-length waveguide 
takes values uniformly distributed between zero and one, 
for nominally identical samples and arbitrarily small vari- 
ations of the frequencyJ^ It is therefore important to 
accurately characterize the radiative coupling mechanism 
between distant QDs in a disordered PHC structure. The 
formalism presented here only relies on the knowledge of 
the spectrum and field maps of the PHC modes. To- 
gether with the possibility to simulate photon modes in 
very long PHC structures - offered by the recently de- 
veloped Bloch-mode expansion method*^ - it therefore 
represents the election method to carry out a systematic 
study of disorder effects on the radiative properties of 
QDs in PHC structures. 
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